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We model polydomain liquid- crystal elastomers by extending the neo-classical soft and semi-soft 
free energies used successfully to describe monodomain samples. We show that there is a significant 
difference between polydomains cross-linked in homogeneous high symmetry states then cooled to 
low symmetry polydomain states and those cross-linked directly in the low symmetry polydomain 
state. For example, elastomers cross-linked in the isotropic state then cooled to a nematic polydo- 
main will, in the ideal limit, be perfectly soft, and with the introduction of non-ideality, will deform 
at very low stress until they are macroscopically aligned. The director patterns observed in them 
will be disordered, characteristic of combinations of random deformations, and not disclination pat- 
terns. We expect these samples to exhibit elasticity significantly softer than monodomain samples. 
Polydomains cross-linked in the nematic polydomain state will be mechanically harder and contain 
characteristic schlieren director patterns. The models we use for polydomain elastomers are spa- 
tially heterogeneous, so rather than solving them exactly we elucidate this behavior by bounding 
the energies using Taylor-like (compatible test strain fields) and Sachs (constant stress) limits ex- 
tended to non-linear elasticity. Good agreement is found with experiments that reveal the supersoft 
response of some polydomains. We also analyze smectic polydomain elastomers and propose that 
polydomain SmC* elastomers cross-linked in the SmA monodomain state are promising candidates 
for low field electrical actuation. 

PACS numbers: 



I. INTRODUCTION 

Liquid crystal elastomers (LCE's) are rubbery net- 
works of cntropically dominated polymer chains that ex- 
hibit mobile liquid crystalline order. The elasticity of 
liquid crystal elastomers is a rich and subtle subject be- 
cause the liquid crystal order is anisotropic, the strains 
involved are often large and the elastomer responds very 
softly to some stresses but more conventionally to oth- 
ers. Monodomain liquid crystal elastomers (elastomers 
in which the liquid crystal director does not vary from 
point to point in the sample) were first successfully syn- 
thesized in 1991 [l|. Since then many remarkable phe- 
nomena have been observed in nematic elastomers - 
extremely soft elasticity and textured deformations Q, 
dramatic spontaneous deformations at the transition be- 
tween the isotropic and aligned states [1| , rotation of the 
nematic director by the application of stress [H or elec- 
tric fields Elongations go back to molecular shape 
depending on the magnitude of nematic order. Mechani- 
cally soft response depends on rotating the director with 
concomitant shears and elongations that accommodate 
the rotating long direction of the chain distribution 0- 

, summarised in 0] . The very special deformations that 
accompany soft response of monodomains depend criti- 
cally on the direction of the initial director to the overall 
extension direction. It therefore seems remarkable that 
polydomain elastomers are also capable of soft response 
to macroscopic deformations [l(J El that also progres- 
sively align the directors. One would expect neighbour- 
ing regions with initially different directors and therefore 
differing soft strains would suffer mechanical incompati- 
bility. Differing strategics for overcoming these apparent 



geometric impediments to soft response have been pro- 
posed [IM3; we discuss this work in detail later. 

In this paper we attack the central geometric issues of 
how extremely soft response is attained in isotropic gene- 
sis polydomain elastomers and why nematic genesis elas- 
tomers are harder. We bound the effect of semi-softness 
in the former response. Recent experiments ficj ] show 
remarkable, genesis-dependent, super-softness of poly- 
domain nematic elastomers. The different response of 
elastomers of differing genesis is a result we share with 
Uchida, though our route to geometric compatibility is 
rather different. These results may be significant to the 
development of LCE actuators since polydomain elas- 
tomers are significantly easier to synthesize than their 
monodomain counterparts, can exert much larger forces 
because they are not limited to the thin film geometries 
needed for making monodomains, and can respond to 
external stimuli, such as polarised light, without having 
such large internal (mechanical) fields to overcome. 



A. Models of Liquid Crystal Elastomers 

Microscopically liquid crystal elastomers can be 
thought of as the result of cross-linking a liquid- 
crystalline polymer melt. The melt consists of long 
writhing polymer molecules that incorporate many short 
rod-like units - fig. [TJ Below a transition temperature 
the rod like mesogens align to form a liquid crystal 
phase. The introduction of cross-links between the poly- 
mer chains prevents them from flowing past each other 
making the material a solid elastomer rather than a liq- 
uid melt. However, sufficiently low density cross-linking 



2 




FIG. 1: Diagram of a nematic polymer melt. The black lines 
represent polymer strands and the ellipsoids are nematic rods. 




FIG. 2: A liquid crystalline polymer is cross-linked in the high 
temperature isotropic state - centre, bottom row On cooling, 
a nematic phase is formed and the elastomer stretches along 
the new nematic director. Since any director could have been 
chosen there are many equivalent low energy states, shown in 
an arc around the isotropic state. Deformations around the 
arc will not cost energy. 

does not stop the extensive thermal motion of the chains 
or prevent the formation of liquid crystal order, so the 
resulting solid is a liquid crystalline elastomer. The novel 
elastic properties of a liquid crystal elastomer can most 
easily be understood by considering an elastomer that 
was cross-linked in the high temperature isotropic state 
then cooled to a nematic monodomain in which the ne- 
matic director is constant throughout the sample. After 
the nematic phase has formed the polymers are more 
likely to adopt conformations which are extended along 
the nematic director, so the elastomer undergoes a uni- 
axial stretch along the nematic director - fig. [2] However, 
since the high temperature cross-linking state is com- 
pletely isotropic, any direction could have been chosen 
for the nematic director so there is a degenerate set of low 
energy states each with a different director and a different 
deformation with respect to the cross-linking state. De- 
formations that take the elastomer from one low energy 
state to another must not cost any energy to impose and 
result in rotation of the liquid crystal director through 
the elastomer. This symmetry-driven phenomenological 
explanation of soft elastic modes was first presented by 
Golubovic and Lubensky [lj|; see [§] for the theorem's 
application to liquid crystal elastomers. 

Experimentally elastomers with spontaneous elonga- 



tions at the isotropic-nematic transition by factors of 4 
[3] have been synthesized, so it is clear that any theory of 
liquid crystal elastomers must be fully non-linear in the 
strain tensor. Moreover, strains induce director rotations 
of up to 90°, hence finite rotation theory is also needed. 
The simplest model that we can write down for the sys- 
tem discussed above is simply a neo-Hookean free energy 
that incorporates a large spontaneous deformation 

F=^iTr{f -T 1 -g) (1-1) 

where [i is a shear modulus, 7 is the deformation gradi- 
ent from the cross-linking state (here the isotropic state) 
and, if <5 is the identity matrix and r is a constant material 
parameter, I is derived from the final state liquid crystal 
director n as I = r~ 1 / 3 (S + (r — l)nn) (a notation differ- 
ing from Q by a factor of r 1 / 3 ). In the isotropic phase 
r = 1 and this form reduces to a standard neo-Hookean 
energy, but on entering the nematic phase r > 1 and 

this energy is minimized by any choice of n provided it 

1/2 

is accompanied by a spontaneous deformation 7 = £ ' 
which is an elongation from the isotropic state along n 
by a factor of r 1 / 3 . Elastomers that suffer spontaneous 
deformations by factors between 1.05 and 4 have been 
observed, implying values for r between 1 and 60 [|[ . As 
for all elastomers, LCE's are essentially incompressible so 

we require that Dot ^7^ = 1, a condition satisfied for ex- 
ample by spontaneous deformations, Det^ 1 ^ 2 ^ = 1. We 

can rewrite this energy in terms of deformations from a 
minimum energy, relaxed nematic state with director no 
by substituting j n = 7 • £q ' to give 

F = i^Tr (jg • ■ I" 1 ' ^) 1 ( L2 ) 

where n is the nematic director in the reference state 
and n is the director in the final state. A microscopic 
derivation of these free energies is given in [ll| which is 
a simple extension of the classical theory of rubber based 
on Gaussian polymer chains. The microscopic model pro- 
vides one key additional insight — in the idealized limit of 
compositionally homogeneous Gaussian chains the local 
energy does not remember the cross-linking state, so the 
above energy for an elastomer cross-linked in the isotropic 
state then cooled to the nematic state is also appropri- 
ate if the elastomer is cross-linked in the nematic state 
with director no. Real chains are never perfectly Gaus- 
sian or completely compositionally homogeneous, so real 
elastomers cross-linked in the nematic state do remem- 
ber the director at cross-linking and their free energy is 
minimized when their director aligns in this direction. 
This effect is modeled by adding a "non-ideal" term to 
the energy which favors alignment of the director along 
no giving a free energy density such as 

F= ±/iTr (2 ' |£ ' 1 T ' I"' + "1 ' (£- n ono) ■ f • nn) , 

(1.3) 
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the most general form at order A [2Cj. The coefficient 
a determines how far the elastomer deviates from ideally 
soft behavior. 

Liquid crystal elastomers can also be prepared in lay- 
ered (smectic) liquid crystal phases. Elastomers have 
been synthesized in both SmA phases in which the di- 
rector is along the layer normal and SmC phases where 
there is a preferred tilt angle 9 between the director and 
the layer normal. Smectic elastomers are modeled by as- 
suming that the layers deform affinely and that there are 
energy penalties for deformations that change the layer 
spacing or cause the director to rotate away from its pre- 
ferred tilt angle. This means that SmA elastomers do 
not have any soft elastic modes since the layers deform 
affinely and rotation of the director away from the layer 
normal costs energy. However, the elasticity of SmA elas- 
tomers is still a complicated subject because the modulus 
for stretching the inter-layer spacing is much higher than 
the modulus for in-plane stretches so their elasticity is 
essentially two dimensional. SmC elastomers in which 
the director makes a fixed angle with the layer normal 
do still have soft modes because the director can rotate 
in a cone around the layer normal without changing the 
inter-layer spacing or deviating from the preferred tilt 
angle. Accordingly we shall predict unusual mechanical 
response for polydomain SmC elastomers. 



B. Polydomain and Monodomain Elastomers 

An LCE is a monodomain if the director is the same 
at every point in the relaxed elastomer and a polydo- 
main if it points in different directions at different points 
in the relaxed elastomer. Synthesizing monodomains is 
difficult because, in order to make the elastomer choose 
the same director at each point, this direction must be 
imprinted on the elastomer at cross-linking, normally by 
cross-linking under uniaxial stress. If the elastomer is 
cross-linked without any such imprinting, in either the 
isotropic or the nematic state, then it forms a polydo- 
main. Polydomains and monodomains are easily dis- 
tinguished because monodomains are highly transpar- 
ent and exhibit large spontaneous deformations at the 
isotropic-nematic transition while polydomains have no 
macroscopic spontaneous deformation, and are opaque 
in the nematic state (because the gradients in the di- 
rector scatter light) but become transparent if they are 
stretched enough to align the director throughout the 
sample (2lj . 

We extend the monodomain energies in the previous 
section to polydomain energies by allowing n, no and 7 to 

become spatially varying fields. We define 7(x) = Vy to 
be the deformation gradient from the cross-linking state. 
This means that, ignoring non-ideal terms, the energy 
function for a polydomain cross-linked in the nematic 
state will be 



where £q is derived from no(x) which is the nematic di- 
rector at cross-linking. In this case the form of this field 
will be the nematic disclination pattern present in the 
nematic-melt before cross-linking. This energy is very 
different to the monodomain energy from which it was 
derived because it has a very significant x dependence in 
no(x). However, the energy function for a polydomain 
cross-linked in the isotropic state is simply 
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which is exactly the same as the corresponding mon- 
odomain energy, cqn. (|l.lj) . There is no intrinsic spa- 
tial variation in this function, which is not surprising 
since the cross-linking state is completely isotropic and 
homogeneous. However, monodomains and polydomains 
cross-linked in the isotropic state are manifestly differ- 
ent. We propose that this is because of the form of the 
non ideal terms that must be added to this energy. In 
the monodomain case the elastomer is cross-linked in the 
presence of a uniaxial stress which imprints a constant 
preferred direction no (in effect a field) on the elastomer 
so that when it cools to the nematic state it adopts the 
same nematic director everywhere, leading us to a mono- 
domain energy 
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In the polydomain case there is no imprinting so there 
is very little to break the isotropy of the cross-linking 
state and introduce any non-ideal terms at all. However, 
weak mechanisms do exist such as cross-linking molecules 
having a rod like character which impose an additional 
direction locally on the network 22(. Consequently, in 
any finite region the cross-linking rods develop a slight 
average orientation fl2l [23[. These mechanisms permit 
the inclusion of a very small non-ideal term with a spa- 
tially varying preferred direction no(x). It is this dis- 
tinction between the large homogeneous non-ideal term 
in the monodomain case and the small spatially varying 
term in the polydomain case that drives the distinction 
between the two systems. 



C. Polydomains - Macroscopically Hard or Soft? 

In this paper we argue that whether polydomain elas- 
tomers are macroscopically hard or soft depends on the 
relative symmetry of the cross-linking state and the final 
polydomain state. In this section we will give a quali- 
tative overview of the cause of this behaviour, then in 
subsequent sections we will work through three examples 
of possible polydomain systems. 

An ideal liquid crystal elastomer energy typically lo- 
cally has a set of deformations that minimize the free- 
energy. These are generated by a symmetry breaking 
phase transition from a high temperature parent state 
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FIG. 3: Left - the ring of deformations leading to low energy 
states for a liquid crystal elastomer. The dot in the center 
is the high symmetry parent state from which the broken- 
symmetry, low energy states on the ring are derived. Right - 
visualization of the quasi-convex hull of the set of low energy 
deformations. These deformations, inside the ring (shaded 
grey), can be imposed with minimal energy because they can 
be constructed out of textures of the soft deformations that 
each lie on different points of the ring. 



that is accompanied by a deformation. Since the tran- 
sition breaks a symmetry there are many ground states, 
each with a different deformation with respect to the par- 
ent state. We visualize this set of energy minimizing de- 
formations as a ring with the parent state at the centre 
— fig. [3] Such a system is vulnerable to the formation of 
textured deformations since, if a deformation on the inte- 
rior of the ring is imposed this is not a low energy state, 
but it is possible that the energy of the deformation can 
be reduced to zero by the elastomer splitting into many 
small regions each of which undergoes one of the low en- 
ergy deformations in such a way that the macroscopic 
deformation is what was imposed. The ability of LCE's 
to form such textured deformations has been the subject 
of several studies — [3, [24T - [26j . We represent the set of 
all deformations that can minimize the energy after the 
formation of the most advantageous textured deforma- 
tions, the quasi-convex hull (QCH) of the set of energy 
minimizing deformations, as the interior of the ring of 
energy minimizing deformations — fig. 

A polydomain sample cross-linked in the high symme- 
try (high temperature) state is effectively cross-linked at 
the centre of the quasi-convex hull and, on cooling to the 
aligned state, forms the same quasi-convex hull at each 
point. Although at each point the elastomer has under- 
gone a deformation that takes it to the boundary of the 
set, these deformations are put together in an elastically 
compatible way so that the elastomer as a whole is in 
a textured state at the centre of the set. Such an elas- 
tomer can be deformed macroscopically softly simply by 
moving the constituent domains around the quasi-convex 
hull, which will cause different soft textures to evolve as 
the elastomer as a whole moves across the quasi-convex 
hull. 

In contrast, if the elastomer is cross-linked in the low 
temperature, symmetry-broken state then, although each 
point still has a quasi-convex hull of the same form each 
domain sits on the edge of its quasi-convex hull and the 
deformation required to take each domain back to the 
center of its hull is different. This corresponds to a pic- 



FIG. 4: If the elastomer is cross-linked in the low symme- 
try polydomain state then, although every domain still has 
a quasi-convex hull of the same form, each domain is cross- 
linked at the boundary of its quasi-convex hull and requires 
a different deformation to get to the its centre. The quasi- 
convex hulls of a few domains are illustrated in the diagram, 
and the cross represents the cross-linking configuration. The 
cross is the only point in every quasi-convex hull so if the elas- 
tomer deforms away from this point some domains will leave 
their convex hull and the response will be hard. 



ture like fig. 2] Although each domain has a hull and 
soft modes, there are no deformations that are soft for 
all domains so the elastomer is macroscopically hard. As 
discussed in the previous section the distinction between 
a polydomain cross-linked in the high symmetry state 
and a monodomain is driven by the addition of differ- 
ent non ideal terms. This is easily visualized in terms of 
the above sets - the ideal monodomain and high symme- 
try cross-linking polydomain have the same quasi-convex 
hull at each point, but the addition of non- ideal terms 
breaks the degeneracy of the states in the quasi-convex 
hull, making one a unique global minimum. In the mon- 
odomain case this is a point on the boundary of the set 
with a uniform director throughout the sample. In the 
polydomain case, the non-ideal terms vary spatially such 
that a different low energy state on the boundary of the 
QCH is favoured for different points. The global mini- 
mum is the point at the centre of the set, no deformation, 
and is achieved by a complicated textured, low energy, 
deformation of the different regions of the polydomain. 

The above analysis of soft polydomains has a very sim- 
ple physical interpretation. If an elastomer is cross-linked 
in the high symmetry state then cooled to the low sym- 
metry state there is no energetic penalty (except small 
deviations from ideality) to stop it breaking the sym- 
metry in the same way at every point in the elastomer. 
Therefore the macroscopic deformations that take the 
polydomains into these well- aligned states must be soft. 
This will not be true if the crosslinking is in the low sym- 
metry polydomain state. 

In the remainder of this paper we will analyze nematic 
polydomains cross-linked in the isotropic state and the 
nematic state and SmC polydomains cross-linked in the 
SmA monodomain state. The first and last of these are 
examples of soft polydomains, while the middle is a hard 
polydomain. 
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II. NEMATIC POLYDOMAINS 

A. Formulating the elasticity problems 

Following the discussion in the introduction, we ex- 
pect a profound difference between nematic polydomains 
cross-linked in the isotropic and nematic state. We 
further distinguish between ideal and non-ideal polydo- 
mains. Ideal systems have no locally preferred director 
orientation and therefore locally completely soft modes. 
The addition of non-ideal terms slightly favours a par- 
ticular local director orientation making the modes that 
were previously completely soft now cost a small amount 
of energy. This leads us to four types of nematic polydo- 
main corresponding to four different free energies, 
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where il/nl denote ideal/no n- ideal elastomers crosslinkcd 
in the isotropic state, and iN/nN ideal/non-ideal elas- 
tomers crosslinked in the nematic state. The bulk mod- 
ulus of elastomers is several orders of magnitude higher 
than the shear modulus, so all deformations are volume 

preserving, that is Dct^7^ = 1. The preferred direction 

no is discussed below. 

Having identified four types of polydomain, ideal and 
non-ideal, with either a nematic or an isotropic crosslink- 
ing state (genesis), we wish to study the energetic cost 
of imposing macroscopically homogeneous stretches on 
large blocks of these different types of polydomains. By 
large, we mean large enough to contain very many do- 
mains and be macroscopically isotropic. In each case 
we take the cross-linking state of the elastomer as the 
reference configuration from which deformations will be 
measured, and define the displacement field from this 
state as y(x), so that the local deformation gradient is 
7(x) = Vy. We then wish to study the energy of a large 

sample occupying a domain f2 (with boundary 5fl) in the 
reference configuration, that is subject to a macroscopic 
deformation A after it has adopted the most favourable 
internal deformation and director pattern. We define this 
relaxed energy function as 

^ r (A,n (x)) = min min — — /F(Vy, n(x), n (x))dx, 

— y(x)s.t. n(x) Vol ll J 

y=A x 

on cSQ 

(2.8) 

and the four different types of polydomain correspond to 
four different choices for F(Vy, n(x), no(x)). 

There is an important distinction between the two 
fields rio(x), which defines the local preferred nematic 
alignment, and n(x) which is the nematic field after de- 



formation. The former is a fixed field for a given poly- 
domain that encodes all its spatial heterogeneity, while 
the latter is a variable field which the elastomer will ad- 
just to minimize its free energy as it evolves under the 
macroscopic A. Furthermore, in the nematic genesis case 
we expect the form of fio(x) to be a disclination texture, 
while in the isotropic genesis case it will be a random 
field arising from the sources of disorder in the crosslink- 
ing state. In the latter case no does not correspond to 
the equilibrium director pattern at zero stress, but rather 
locally it is the director that a domain would adopt if it 
were unconstrained by its neighbors (e.g. by cutting it 
out of the sample). These locally optimal strain fields 
7(x) associated with this director pattern are extremely 
unlikely to be compatible deformations (be the gradient 
of a continuous displacement field) and thus would re- 
quire the sample to fracture. 

Our task is to find the lowest energy, compatible strain 
fields - in general a difficult task. It is not clear that 
they are very subtle fields associated with a continuous 
director variation (but still compatible) or are compati- 
ble combinations (textures and laminates) of strains as 
are observed in the deformation of monodomains under 
boundary constraints at variance with soft deformations 
[H . We shall proceed to find bounds on the elastic energy. 
In the ideal case, we know that textures of deformations 
each on the boundary of the QCH can give zero energy 
cost for macroscopic deformations within the QCH and 
thus give an exact value for the energy (a continuous 
field could not do better than this choice of test field). 
In the non-ideal case we will use these exact minimizers 
of the ideal part of the free energy and evaluate the ex- 
tra, non-ideal cost associated with them, thus forming an 
upper bound. In reality an elastomer could, by adjusting 
the laminates, or by finding an unlikely continuous field, 
lower the energy by correlating distortion fields with the 
random field. We discuss that energy reduction strategy 
briefly later. 



B. Ideal isotropic genesis polydomains 

The only completely solvable system is that of ideal 
isotropic genesis (il). The energy function has no no de- 
pendence so it is not really spatially heterogenous at 
all. This is unsurprising since the isotropic cross-linking 
state appears to be completely homogeneous. It is the 
same model as is used for ideal monodomain samples, 
so it has been intensively studied. Its key property is 
that it is minimized locally by deformations of the form 
7 = Vy = i 1 ^ 2 , uniaxial stretches by r 1 / 3 along n which 
can be in any direction. These correspond to the spon- 
taneous deformations that are caused by the transition 
to an aligned state. Since there are no preferred direc- 
tions for alignment, one possible deformation field that 
allows every domain to be stretched uniaxially by r 1 / 3 is 
simply for the whole sample to stretch this much along 
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the same axis, so the deformation gradient is Vy = I 1 
at every point in the sample. This is a macroscopically 
aligned, monodomain, state with the director along the 
stretch axis and is a minimum energy state on the bound- 
ary of the quasi-convex hull. Clearly any direction could 
have been chosen for the stretch, so there is a degenerate 
set of equivalent low energy states each with a different 
deformation with respect to the cross-linking state, and 
deformations that map between these states must be soft. 

Other macroscopic deformations will also be minimum 
energy states if they can be realized as the average of a 
compatible deformation fields that are uniaxial, volume- 



preserving stretches by 



,1/3 



at every point in the elas- 



tomer - these macroscopic deformations lie on the inte- 
rior of the quasi-convex hull. However, the condition of 
mechanical compatibility places a very strict condition 
on these deformation patterns. If two adjacent regions 
undergo different uniaxial stretches then this will gener- 
ally lead to the material fracturing along the boundary 
between the regions. However, if the regions are sepa- 
rated by a plane boundary and the boundary bisects the 
angle between the two stretches, then the boundary will 
be stretched to the same degree by both deformations 
and material continuity can be restored by simply body 
rotating both regions. The fracturing and compatible 
possibilities are illustrated in Fig. [5] This latter type of 
deformation pattern can be repeated many times to make 
a textured deformation that achieves a total deformation 
that is not simply a uniaxial stretch by r 1 / 3 whilst ensur- 
ing that the deformation at every point is such a stretch. 
The structured type of director and deformation pattern 
this leads to - oscillations between regions of constant de- 
formation and director separated by plane boundaries - 
is characteristic of textured deformations that arc driven 
by the requirement of mechanical compatibility between 
domains [3, [I?} . The stripe-domains observed in mon- 
odomain samples (28l . [29j are another example of this 
type of behavior. 




FIG. 5: Top: Two uniaxial stretches with different axes ap- 
plied either side of a boundary generally result in the material 
fracturing at the boundary. Bottom: If the boundary normal 
bisects the two axes then the boundary is stretched to the 
same degree under both deformations and material continu- 
ity can be restored by body rotating the two regions back 
together. 

The full relaxation, F r (X), of the ideal isotropic gen- 
esis energy has already been found [24j in the context 
of ideal monodomain elastomers. Remarkably, there are 



textured deformations that allow any macroscopic de- 
formation with largest principal stretch less than r 1 / 3 
and smallest principal stretch greater than r -1 / 6 to be 
achieved whilst every point in the material undergoes a 
local uniaxial stretch of r 1 / 3 . This class includes the de- 
formation A = 5, so the elastomer can accommodate its 
spontaneous deformation without changing shape macro- 
scopically, and every uniaxial stretch up to r 1 / 3 , so the 
elastomer can be stretched by any degree up to r 1 / 3 with- 
out the energy of the elastomer rising. After a uniaxial 
stretch of r 1 / 3 the elastomer is in a completely aligned 
monodomain state, so further stretching in the same di- 
rection cannot lead to further director rotation and the 
elastomer deforms without texture and is no longer soft. 

The full relaxation result in [24[ is that, if the principal 
stretches of A are /i < /2 < /sj 
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K qc = {A G M 3x3 : r~ 1/6 < fi < fa < r 1/3 , Dct(A_) = 1}. 

The "else" case only contains deformations with 
Oct (A) 7^ 1 that do not conserve volume. Textured de- 
formations are required for macroscopic deformations in 
/ and K qc but not S. The set K qc is an eight dimensional 
set of zero energy deformations, so if an elastomer is at 
a point in the interior of K qc all small volume preserving 
deformations of the elastomer are also in K qc and do not 
cost energy to impose - the the energy is liquid-like. In S 
the energy depends on all three principal values of the de- 
formation so the energy is solid like while in / it depends 
on only the smallest principal value so the energy is in- 
termediate between that of a solid and that of a liquid. 
Physically, the set S consists of stretches, fa, sufficient 
that the director has completely aligned with this stretch 
direction, so the elastomer responds to further stretching 
in the same way as a conventional rubber, while the set 
K qc consists of deformations that can be made soft by the 
formation of textured deformations of spontaneous defor- 
mations. The set / contains compression in one direction 
making the sample thinner than any spontaneous defor- 
mation would. The director lies perpendicular to the 
compression axis but still forms textures in that plane 
allowing any macroscopic in-plane shapes bounded by 
fal fi < \fr. The energy however only depends on the 
degree of compression. 

In this work we are concerned with uniaxially stretch- 
ing polydomains to induce the polydomain-monodomain 
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transition. This corresponds to macroscopic deforma- 
tions of the form A = diag(A, l/v% 1/vA) for A > 1. 
Applying the above relaxation result, these deformations 
will be achieved at the energies 



Fn(X) 



3/£ 



if 1 < A < r 1 / 3 
if A>r!/3. 



(2.10) 



Differentiating this gives the engineering stress, 



0"ir(A) 



dFn 



dX 



JO A < rV3 

[fir^iX/r- 1/A 2 ) A>r!/3. 

(2.11) 

This simply means that ideal isotropic genesis polydo- 
mains will deform at zero stress but with textured defor- 
mations until they have been stretched by r 1 / 3 , at which 
point they are completely aligned monodomains and re- 
spond to further deformation as a neo-Hookean solid. 
This result is quantitatively wrong, real isotropic gen- 
esis polydomains do not deform a absolutely zero stress, 
and when stress is removed they (macroscopically) re- 
turn to their original configuration. This motivates the 
consideration of non-ideal theories. However, since the 
cross-linking state is almost isotropic, non-ideality must 
be very small so we expect several key features of the 
ideal model to persist. Extension will not occur at zero 
stress, but extensions up to A < r 1 / 3 will occur at ener- 
gies 0(a) <C 1. Furthermore, the deformation patterns 
and director patterns in the non-ideal case must still be 
very close in energy to those in the ideal case, so the 
observed patterns will still be characterized as textured 
deformations driven by clastic compatibility - not discli- 
nation textures. Indeed, nematic disclinations never have 
zero-energy elastically compatible associated distortions 
[30j . so cannot be observed in isotropic genesis systems. 



C. Non-Ideal isotropic genesis polydomains 

Finding the full relaxation of the non- ideal polydomain 
energy is probably intractable. At the moment the relax- 
ation is not known for the easier monodomain case (ex- 
cept in a thin film limit [25j ) , and the polydomain result 
will depend to some extent on the exact form of no(x). 
However, we can put upper and lower bounds on the 
energy-strain curves. Developing an upper bound on the 
energy is straightforward. We simply use a textured test 
strain field from the ideal case and calculate its energy 
in the non ideal case. Since the relaxed energy function 
is a minimum over all strain fields, evaluating the energy 
at one example of a strain field is an upper bound on the 
energy. In a sense this is a Taylor-like bound of uniform 
strain, but our uniform macroscopic strain is in fact com- 
posed of textures that allow deformation anywhere in the 
QCH. We depart also from conventional Taylor bounds 
in that our bound is valid for large strains, up to r 1 / 3 
which can be 100s% for nematic elastomers. 



If a point in the elastomer undergoes a uniaxial exten- 
sion of magnitude 7 from the crosslinking state, at an 
angle 9 to the preferred direction no, then the energy of 
the deformation is 



F = min 



1111 — Ir 

n 2 V: 



7 -7 T \ 



i/xr 1 / 3 



1 + ar 1 / 3 7 (6_ — n n ) j T iinj 
(2.12) 



where the minimization over n is achieved by taking n 
along the axis of 7. If the region of constant deformation 

is much larger than the individual domain size (region of 
given no) then averaging over no gives: 



(2.13) 



which is minimized at 7 3 n = r/(l + 2ar/3) with a value 



W 2_ar 
2 V 3 



1/3 



(2.14) 



This free energy density can be achieved if every point 
in the sample undergoes a uniaxial elongation of magni- 
tude 7 m in any direction. This situation is completely 
analogous to the situation in ideal elastomers in the 
isotropic configuration where the same (in this case min- 
imal) energy can be reached by applying an elongation 
of magnitude r 1 / 3 in any direction. The DeSimone and 
Dolzmann texture result (2~H shows that any uniaxial 
macroscopic deformation with magnitude less than 7 m 
can be achieved by a texture of deformations in which 
each deformation is a uniaxial deformation by j m . This 
allows us to place an upper bound on the total energy of 
the sample after it has undergone a macroscopic uniaxial 
elongation by A, 



2F r (X) 1 3(1 



2ar \ 1 / 3 
3 



A 3 < r/(l+2w/3) 



H jrV3 (2 + (I + 2a) A 2 ) A 3 >r/(l + 2ar/3). 

(2.15) 

Although this upper bound has been calculated by using 
textures with regions of constant deformation that are 
very large compared to the length scales no varies on, 
the same result can be achieved with any size. This is 
because the averaging to 1 /3 of the cos 2 9 in cq. (|2.12[) 
will still be true after averaging over many domains pro- 
vided the axes of the domains are not correlated with the 
no field. Introducing such correlations would reduce the 
energy of the elastomer and would also determine length 
scales involved in the actual deformation and director 
fields. An Imry-Ma style attack on this problem, but not 
involving textured test fields as here, is due to Terentjev 
and Fridrikh [l2[ . Domain size would be selected to take 
advantage of fluctuations in the random ordering field 
from crosslinking. We return to this problem elsewhere. 

The ideal system provides a very simple lower bound 
on the energy, eq. (|2.10[) . since the non-ideal term is 



never negative. We can improve on this bound by us- 
ing the Sachs limit of stress uniform through the sam- 
ple. This provides a lower bound because it neglects 
the requirement of compatibility of deformations (7 = 

Vy). We calculate this bound numerically by minimiz- 
ing F„j yy, n, ht^j — aj xx across all 7 and n at fixed a for 
a given domain (no) to find the optimal deformation 7 

=m 

and director orientation n m of the domain at the stress 
a. The energy and extension of the whole sample are 
then found by averaging F n j I 7 , n m , no ) and ( 7 ) 

V=m / \=mJ xx 

across all domain orientations. 

Although wc have calculated the full Sachs bound nu- 
merically, we can understand its behavior at small exten- 
sion analytically. At zero stress every domain is free to 

undergo an energy minimizing spontaneous deformation 

1/2 

- anything of the form 7 = R ■ £ ' , where R is a rota- 

• 1/2 
tion. The xx component of this is simply x • R ■ £ ' • x, 

which is to say it is the component of 7 ■ x that is par- 
allel to x. For each domain a rotation R can be chosen 
such that A lies anywhere between and {lo 1 ^ 2 ■ x|. 
Since we are studying extension we are not concerned 
with X xx < 1, but this means that any X xx between 1 and 
(^(q 1 ^ 2 • can, in the Sach's limit, be achieved without 
stress. Therefore the Sach's limit on the energy is simply 
3/i/2 for 1 < X xx < ({io 1 ^ 2 ■ x|Y It is straightforward to 
calculate this average giving 



3.5 



1/2 



„-l/6 



sinh 1 \/r — 1 



sfr 



1 



(2.16) 



This threshold tends to l+0(e 2 ) for r = 1+e, so for small 
anisotropics it is insignificant, but for large r it tends to 
^3/2. 

The energy plot, Fig. El shows that the bounds con- 
strain the energy very tightly: the bounds reproduce the 
exact ideal result if a = 0, and a is expected to be very 
small for isotropic genesis. The Sachs free energy is in 
effect plotted parametrically since one sets the stress and 
obtains the free energy and strain by minimisation and 
averaging. The bounds show that stretching the elas- 
tomer by ~ r 1 / 3 cannot require an increase in the energy 
density of more than 3^/2((l + 2ar/3) 1 ' 3 — 1) ps fiar/3. 
This means that although the extensions up to r 1 / 3 can 
take place at finite stress, the stress cannot be higher than 
dF/dX ~ nar/mr 1 ' 3 - 1)], which will be a very small 
number since a is small. As in the ideal case, extensions 
larger than ~ r 1 / 3 behave in a nco-Hookean manner. 

Wc can calculate bounds on the stress-strain curve for 
these extensions by using the requirement, which ap- 
plies to all one dimensional elastic energies, that the re- 
laxed energy curve be convex, meaning the stress curve is 
monotonic for A > 1. This means that at a given exten- 
sion X xx not only must the energy function lie between 
the two bounds, but the gradient of the energy function 
must not be so great that, if extrapolated forwards as 



2F/ju 



,3(l+2ccr/3) 



.1/3 




1 



,1/3 



FIG. 6: Bounds on the free energy density of an isotropic 
genesis polydomain nematic elastomer as a function of strain, 
plotted in units of ^r 1//3 /2 with r = 8 and a = 0.005. Upper 
curve (dashed, blue): upper bound from a test strain field. 
Middle curve (smooth, red): Sachs lower bound. Lower curve 
(dashed, green): ideal result, lowest bound. 




FIG. 7: Upper and lower bounds on the stress-strain curve for 
a non-ideal isotropic genesis polydomain. Material Parame- 
ters r — 1.65 and 01 = 0.01 and /1 = 33000. The middle curve 
is the derivative of the Sachs lower bound on the energy, which 
provides as estimate of the stress-strain curve, while the up- 
per and lower curves are the bounds derived by requiring the 
energy function be convex. The circles are stress-strain data 
for a real sample. We thank K. Urayama for permission to 
reproduce this data. 



a straight line, the energy curve intersects the upper- 
bound (used above for estimating the approximate max- 
imal stress), or so little that if extrapolated backwards it 
does the same. This bounds the gradient of the energy 
and hence the stress at each extension. 

A plot of the two bounds on the stress strain curve 
is shown in Fig. [3 which shows that the very soft stress 
plateau for A < r 1 / 3 does indeed survive the introduction 
of semi-softness. The experimental results of Urayama et 
al for isotropic genesis polydomain stress [Io| are also 
shown and display pronounced softness - clearly poly- 
domains can deform softly and the requirement of com- 
patibility between domains does not appreciably harden 
response. 

Random field models of polydomain nematic clas- 
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tomers [12Hl6j address the question of how polydomain 
elastomers form, and how the domain structure varies 
with temperature and stress. Uchida, in 2-D numeri- 
cal simulations of domains, found disclinated structures. 
He also concluded that nematic elastomers cross-linked 
in the isotropic state should display soft-behavior while 
those cross- linked in the nematic state would not, using 
a combination of infinitesimal analysis and finite strain 
simulation, both also in 2-D. Our finite-strain analysis 
has been very different: we employ multi-scale test tex- 
tures giving identically zero upper and lower bounds on 
elastic energies (for ideal isotropic genesis elastomers) 
and low energy response in the presence of disorder. In 
common with Uchida, we obtain hard response for ne- 
matic genesis elastomers: 

D. Nematic genesis polydomains 

We expect the coefficient of non-ideality, a, to be 
much higher for nematic genesis polydomains because 
the cross-linking state is not isotropic and distinguishes 
the direction of the nematic director, no. It will proba- 
bly have a similar magnitude to that observed in mon- 
odomain elastomers where a ~ 0.1. However, although 
non-ideality will be larger, it is conceptually less im- 
portant because already the ideal part of the energy 
will contribute significantly to the stress. This is be- 
cause even an ideal nematic genesis polydomain is not 
expected to have any soft modes. Soft modes are gen- 
erated by symmetry breaking spontaneous distortions 
from an isotropic reference state. Locally an individual 
domain does have an isotropic reference configuration, 
which is reached by applying the inverse spontaneous de- 
formation 7 = £q ~ 1 ^ 2 , but the spatially-dependent de- 
formation 7(x) = £o _1 ^ 2 ( x ) w iH n ot be mechanically 
compatible, so it is impossible to apply a deformation 
that places the whole sample in the isotropic reference 
state simultaneously. The reason for this large differ- 
ence between nematic and isotropic genesis polydomains 
is that there are very few deformation patterns of the 
form A(x) = £o _1//2 ( x ) that are compatible deformations, 
and no nematic disclination patterns are compatible de- 
formations. An isotropic genesis polydomain is forced 
to undergo such a deformation on cooling to the ne- 
matic after cross-linking, so it must choose one of the de- 
formation patterns that is mechanically compatible and 
these, as we have seen in the previous section, allow for 
macroscopic movement across the QCH. In contrast the 
nematic-genesis polydomain undergoes its "spontaneous- 
distortion" in the melt where there are no conditions of 
compatibility because adjacent regions can flow past each 
other, a deformation which would result in a fracture in a 
solid network. One consequence of this is that when a ne- 
matic genesis polydomain is heated to what would be its 
"isotropic" state, it is unable to undergo its energy mini- 
mizing contraction 7(x) = (x) everywhere because 



this is not a compatible deformation. This will result in 
the high temperature state being internally stressed, and 
may lead to elevation of the transition temperature. 

We can bound the ideal nematic genesis free energy 
in the same way we bounded the isotropic genesis non- 
ideal free energy. An upper bound is provided by a test 
strain field. The simplest test field would be a constant 
strain throughout the sample, giving a Taylor bound. 
However, we can find a tighter upper bound by apply- 
ing the same strain to each domain, but allowing each 
domain to form textured deformations that average to 
the strain imposed on it; this is what we termed in the 
previous section a "Taylor-like bound". Consequently, 
if the strain imposed macroscopically is A, the free en- 
ergy density of each domain will be -Ff^A • i 1 ^ 2 ), where 
Ffi(X) is the relaxed energy for an ideal monodomain in 
[24] given by eq. (|2.9|) and is attained by texturing (by 
laminates). The factor of l 1 ^ 2 in the argument of 
is appropriate because the function Fj^X) is written in 
terms of deformations from the isotropic reference state, 
whereas the domains in the nematic genesis polydomain 
are already in the elongated nematic state. The defor- 
mation £ 1 is the deformation the isotropic state would 
have to undergo to reach the nematic state that the do- 
main was cross-linked in, and the component A, of the 
compound deformation in the argument of F^ , is the de- 
formation from this nematic state. A lower bound can 
be found using the Sachs constant stress limit. However, 
since the elastomer is ideal, the individual domains can 
deform softly. This means that, in the constant stress 
limit where compatibility of the deformations between 
the different domains is not required, the elastomer can 
deform completely softly. We can calculate the end of 
the soft plateau in the Sachs bound analytically in the 
same way we did for the isotropic-genesis case. Since one 

1 /2 • 

considers a spontaneous contraction £ in going from 

the nematic state to an isotropic reference state, followed 
by a spontaneous elongation r 1 / 3 in elongating to the ne- 
matic state along n, the end softness in the Sachs bound 

will occur at (^r 1 / 3 ^ which evaluates to 

This limits to 1 + e/4 for small anisotropy r = 1 + e, and 
to y/rir/4 for large r. 

The two bounds on the energy (calculated numerically) 
are both plotted in Fig. [9] These bounds on the en- 
ergy are not very good — there is a large gap between 
them. The Sach's limit displays complete ordering of 
the elastomer at zero stress while the test strain field 
shows hard elasticity and finite modulus (~ /i) at all ex- 
tensions. This is because nematic genesis polydomains 
are strongly heterogeneous materials, and good methods 
for finding stress strain curves for such materials have 
yet to be developed at large deformations that are re- 
quired here. Taylor (affine) and Sach's (constant stress) 



FIG. 8: An example of a possible director pattern at cross- 
linking in a nematic genesis polydomain that certainly doesn't 
have any macroscopic soft modes. This can be seen by con- 
sidering a line element along the boundary in the plane of 
the diagram - all soft modes of the vertical stripes require 
this element to contract while all soft modes of the horizon- 
tal stripes require it to extend or not deform, so there are no 
macroscopic soft deformations. 

bounds on heterogeneous materials are independent of 
the domain structure of the material, meaning the Tay- 
lor bound gives an upper bound on the energy of the 
hardest possible domain structure, and Sach's gives a 
lower bound on the energy for the softest possible do- 
main structure. Unfortunately there exist vanishingly 
unlikely domain structures which are indeed completely 
soft, namely the textured-deformation domain structures 
realized by the isotropic genesis polydomains. There also 
exist domain structures which certainly have no macro- 
scopic soft modes - fig. [5J Since the domain structures 
in nematic polydomains will certainly not be of the com- 
patible double laminate type — they will be disclination 
textures which are not mechanically compatible — the 
lower bound is of little physical significance. Therefore, 
we estimate the stress strain relation for these polydo- 
mains by simply taking the derivative of the upper bound 
on the energy. The resulting stress-strain curve, shown in 
Fig. I10[ matches the experimentally observed completely 
hard behavior. 

The inclusion of non-ideality in the nematic genesis 
model will elevate the Sach's limit to a plateau of height 
cx a but will not alleviate the fundamental difference be- 
tween the soft elasticity seen in the Sachs and the hard 
elasticity seen in the Taylor-like limit. 

III. SMECTIC POLYDOMAINS 

Smcctic liquid crystal phases are phases in which the 
rods not only have orientational order but are also lay- 
ered. Liquid crystal elastomers can exhibit smectic or- 
dering [3l| . and monodomains with both SmA ordering 
[3^ | (in which the liquid crystal director is parallel to the 
layer normal) and SmC ordering [33| (in which the liquid 
crystal director makes a constant angle 9 with the layer 
normal) have been synthesized. SmC phases in which 
all the rods have the same chirality (SmC* phases) are 
particularly interesting because they have electrical po- 
larizations along the cross product of the layer normal 
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FIG. 9: Bounds on the free energy density of an ideal nematic 
genesis polydomain with r = 8. Upper curve (green) — upper 
bound using a test strain field. Lower curve (blue) — lower 
bound using a constant stress. There is a very large gap 
between the two bounds. Note that the ideal system displays 
hard elastic response for purely geometrical reasons stemming 
from compatibility requirements. 




FIG. 10: Estimate of stress vs strain for an ideal nematic 
genesis polydomain elastomer with r = 1.65 and [i = 37000, 
obtained by differentiating the upper bound on the energy. 
The circles are stress-strain data for a real sample, the ex- 
act analogue of the sample used in fig. [7] but cross-linked in 
the nematic polydomain state. We thank K. Urayama for 
permission to reproduce this data. We note that the fitting 
parameter r is the same but a slightly higher value of /1 is 
needed. This probably reflects the fact that in reality ne- 
matic genesis polydomains have a significant non-ideal term 
(a ~ 0.1) which will slightly harden the system. 

and the director. The introduction of these phases sig- 
nificantly increases the total number of polydomains that 
can be considered since there are now four distinct states 
— isotropic, nematic, SmA and SmC — and polydomains 
can be made that have been cross-linked in any one of 
these states then cooled or heated to any other of the 
states. 

In strongly coupled elastomers, smectic behavior is 
usually modeled by assuming that the layers deform 



11 



affinely under deformations (so that a layer normal k be- 
comes 7 _T -k after a deformation 7 has been applied) and 

then adding terms to the underlying nematic free energy 
that penalize changing the inter-layer spacing and rotat- 
ing the director away from its preferred angle with the 
layer normal. This means that although SmA elastomers 
have the same cylindrical symmetry as nematics they do 
not have any soft elastic modes since the layers deform 
affinely (they cannot translate or rotate relative to the 
rubber matrix) and rotation of the director away from the 
layer normal costs energy. SmC elastomers could still ex- 
hibit soft modes since their director can rotate in a cone 
around the layer normal without changing the inter-layer 
spacing or deviating from the preferred tilt angle. 

(i) Isotropic Genesis: The modeling assumptions out- 
lined above are well established for monodomain smec- 
tic elastomers but they seem rather strange for polydo- 
main elastomers crosslinkcd in the isotropic state. In 
particular, the assumption that the layers move affinely 
seems appropriate if the layers have been embedded into 
the elastomer at crosslinking. But if, for isotropic gen- 
esis, they have appeared after crosslinking as the result 
of a symmetry-breaking isotropic-SmA transition, then 
they could equally well have formed in any other di- 
rection so one would expect them to be able to rotate 
through the sample. Indeed, since the isotropic-SmA 
transition will be accompanied by a spontaneous defor- 
mation that is a stretch along the director and hence also 
the layer normal, a deformation that returns the elas- 
tomer to its isotropic configuration and then stretches 
it by the same amount in a different direction must, on 
symmetry grounds, be soft and must cause the layers to 
rotate. With this view, the isotropic cross-linked SmA 
polydomains are no different to the isotropic cross-linked 
nematic polydomains since the spontaneous deformations 
at the transitions are of exactly the same form and break 
the same symmetry. Accordingly we expect exactly the 
same macroscopic soft elastic response. SmC polydo- 
mains cross-linked in the isotropic state will also behave 
in the same way. However, although this simple symme- 
try argument cannot be circumvented in equilibrium, it is 
possible that non-equilibrium kinetic effects prevent the 
elastomer from realizing this soft elasticity on experimen- 
tally accessible time scales and that the layers, though 
really free to rotate through the sample, are effectively 
frozen in at the transition to the layered state. This 
freezing would result in a complete hardening of SmA 
elastomers since, while the layers are deforming affinely, 
there are no local soft modes. Equally, affine layer defor- 
mations imposed by the freezing-in of layers would make 
the SmC polydomain much like the nematic genesis ne- 
matic polydomain case since it would possess local soft 
modes but there would be no compatibility between the 
soft modes of adjacent domains so these cannot be used 
to make macroscopic soft modes. 

(ii) SmA or SmC Genesis: Most of the other routes to 
smectic polydomains will result in macroscopically hard 
elasticity, for example cross-linking directly in a SmA 



state will lead to a polydomain with no local soft modes, 
cross-linking in a SmC polydomain state will result in a 
system analogous to the nematic cross linking nematic 
polydomain. 

A. Soft polydomain smectic elasticity 

To recover macroscopically soft elasticity, we need the 
elastomer to break a symmetry at a transition after cross- 
linking in such a way that if it had happened to break it in 
the same way in every domain a monodomain would have 
formed. One interesting system that has this property is 
a SmA monodomain that is cooled without any external 
influences into a SmC polydomain. Crosslinking in the 
SmA monodomain state guarantees that the smectic lay- 
ers are permanently embedded in the elastomer and will 
subsequently deform affinely and that, since the layers 
are embedded as a monodomain, after cooling the system 
can access fully aligned SmC monodomain states without 
having to rotate the layers through the elastomer. To dis- 
cuss this case will take the SmA cross-linking state as the 
reference state and use a set of axes such that the SmA 
layer normal k = (0,0,1) while (0,1,0) and (1,0,0) are 
perpendicular vectors in the SmA layer plane, see Fig. 
[TT1 The local spontaneous deformations at the SmA- 




FIG. 11: A SmA elastomer (left) with director n aligned along 
the layer normal k cools and undergoes a deformation 70 to 

form a SmC elastomer (right) in which the director forms 
an angle 9 with the layer normal. The deformation typically 
includes a contraction along k since the rods have tilted which 
reduces the inter-layer spacing. 

SmC transition if the director tilts in the k — (1,0,0) 
plane will be 

(7n 713 \ 
1/7H733 (3.18) 
733 / 

where all the components of this deformation have fixed 
values determined by the microscopic details of the elas- 
tomer. This transition is illustrated in Fig. [TTJ Elas- 
tomers have been synthesized with values of 713 as large 
as 0.4 H|. 

In the SmA phase there is nothing to distinguish any 
direction in the (0, 1,0) — (1, 0, 0) plane, so the above de- 
formation would also have been soft if it had been applied 
at any other angle in the (0, 1,0) — (1, 0, 0) plane, so the 
full set of soft deformations, Kc; mC , is all deformations 
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that can be written in the form R ■ 70 • i? k , where R is 

a rotation and i?, is a rotation about k. The full set 
— Y 

of deformations that can be made soft by constructing 
textures out of such deformations is known to be [26j 

Kf mC ={A G M 3X3 : Dct(A) = 1, 1 A • k| 2 < 7l3 + 7 | 3 , 
|r T -k| 2 <l/^3,/i(A-^)>l}. (3-19) 

where the function /1 returns the smallest princi- 
pal value of its argument and the deformation A. = 
diag(7i 1 733 , 711733 , p) where 

P= a I'f^k 4 ■ (3-20) 

7l3 + 753 - 713733 

The set Kg C mC is fairly complicated, but the four condi- 
tions in it each have a simple interpretation. The deter- 
minant condition Det (A) = 1 requires that the elastomer 
have the same volume after deformation as it did in the 
SmA state. A line element i n the k di rection of the SmA 
is stretched by a factor of yf 7 2 3 + 733 by any of the local 
spontaneous deformations in K®. mCl so the longest a line 
element in the k direction can be after deformation is 
s/ 7 2 3 + 7 | 3 (which occurs when every local spontaneous 
deformation is the same), imposing the first inequality. 
The second inequality can be understood in an analogous 
way — a vector area along k has its area increased by a 
factor of 1 /733 by any of the local spontaneous deforma- 
tions but its direction may be rotated. The average of 
all these vector areas, A _T • k, cannot be larger than the 

sum of the areas that make it up so |A~ T ■ k| 2 < I/733. 
The final inequality is the analog of /i(A) > 7-" 1 / 6 in the 
ideal nematic case, requiring that the elastomer cannot 
be compressed so much in any direction that it is thinner 
than the natural width of the thinnest direction of the 
underlying chain distribution. 

The set K C g mC is much richer than its ideal nematic 
counterpart, in particular uniaxial deformations of the 
form A = diag(l/v / A, l/-\/A, A), which are stretches by A 
along the original layer normal, are in Kg mC provided 
that 7I3 < A 2 < 733 + 7i 3 • This means that there is a 
whole set of textured SmC polydomain states whose de- 
formation with respect to the parent SmA state is a sim- 
ple stretch along the SmA layer normal. All of these SmC 
polydomain states have the same macroscopic cylindrical 
symmetry as the SmA state. When the SmA sample is 
cooled to the SmC state it could form any of these states, 
so the resting configuration of a SmC polydomain may 
have a uniaxial stretch relative to the SmA state. This 
was not the case in the isotropic genesis nematic poly- 
domains because there was only one textured state with 
the full isotropic symmetry of the cross-linking state. 

In the ideal SmC case our inability to uniquely identify 
one polydomain state with cylindrical symmetry makes 
no difference to the analysis at all since every state in 
KqLq is an energy minimizing state and deformations 
that move the polydomain between them are perfectly 



soft. The addition of a non ideal term, which must be 
small since it breaks both the homogeneity and symmetry 
of the cross- linking state by energetically favoring a single 
director orientation, will have a very similar effect to its 
addition in the nematic case. Non-ideality will break the 
complete energy degeneracy of Kg c placing some states 
(those near the boundary where these is less freedom to 
choose between different textures to minimize the non- 
ideal energy) slightly higher in energy so that small but 
finite stresses are needed to allow the elastomer to explore 
the complete set. The lack of a single unique state with 
cylindrical symmetry means that we can not be sure what 
the energy minimizing state is, indeed which state it is 
will depend on the precise functional and spatial form 
of the non-ideal term included, so it may depend on the 
chemical nature of the elastomer. 



B. Electromechanical switching of soft polydomain 
SmC* elastomers 

Interest in SmC elastomers is mostly driven by their 
potential for electrical actuation. Chiral SmC* liquids 
exhibit an (improper) ferro-electric polarization along the 
cross product of their director and layer normal (35j . Be- 
ing chiral there is a twist of the tilt direction (and hence 
of the polarization) about k on advancing along the layer 
normal direction. There is accordingly no macroscopic 
electrical polarisation unless the twist is undone by ex- 
ternal fields or by boundary effects [36| . On crosslinking 
in the SmA state and cooling to the SmC* state with 
domains, this twist is largely suppressed. Twist of n 
about k means the 713 spontaneous shear direction in 
cqn. (|3.18|) rotates from layer to layer, giving rise to elas- 
tic incompatibility. The effect of such incompatibility 
on texture formation has been discussed in detail in con- 
nection with mechanical switching of SmC monodomains 
[37l ]. Without twist, domains accordingly develop a net 
polarisation so that, when an electric field is applied, en- 
ergy is minimized by domain reorientation so that polar- 
ization is parallel to the applied field. Our inability to 
specify which point deep in the interior of K q s c mC is the 
lowest energy resting state of SmA monodomain genesis 
SmC* polydomains does not prevent us from analyzing 
their electrical actuation since the most extreme actu- 
ation is achieved by making the elastomer traverse the 
whole set Kg c (from boundary to boundary). If an 
electrical field is applied in the (0, 1, 0) direction to such 
a SmC polydomain then this will cause it to form a state 
with its polarization vector uniformly in the (0, 1,0) di- 
rection which it can do by forming a monodomain with 
its director and layer normal both in the (1, 0, 0) — (0, 0, 1) 
plane. The deformation of this monodomain with respect 
to the parent SmA state, 7 , is of the form of eqn (|3.18l) . 

If the electric field is then reversed the elastomer will 
flip into the opposite state which still has the director 
and layer normal both in the (1,0, 0) — (0, 0, 1) plane but 
with the director on the other side of the layer normal 
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so that their cross-product (and hence the polarization) 
is reversed. This state has a deformation with respect to 
the SmA 

(7n ~7i3 \ 
I/711733 . (3.21) 
733 / 

The full deformation undergone be the elastomer when 
the electric field is reversed is A = 7^ • 7 _1 , giving 

/ 1 -2713/733 \ 
|= 1 , (3.22) 

\o 1 J 

simply a reversal of the spontaneous shear. Since the dis- 
crepancy in energy between different states in Kg C mC is 
driven entirely by the addition of non-ideality which ex- 
pected to be small, the electric fields required to perform 
this very large actuation will also be small — smaller than 
the fields used to perform similar actuations on SmC* 
monodomain samples which have large non-ideal fields 
cross-linked into them to make them form monodomains. 
This suggests that SmA monodomain genesis SmC* poly- 
domains are probably better candidates for electrical ac- 
tuation than their monodomain counterparts. 

IV. CONCLUSIONS 

There is a fundamental difference between those poly- 
domain liquid crystal elastomers cross-linked in a high 
symmetry state then cooled to a low symmetry state 
and those crosslinkcd directly in the low symmetry state. 
The former will be extremely soft macroscopically while 
the latter will be mechanically hard. We have ana- 
lyzed two completely soft examples, nematic polydo- 
mains cross-linked in the isotropic state and SmC poly- 
domains cross-linked in the SmA monodomain state and 
one hard example - nematic polydomains cross-linked in 



the nematic state. This distinction between soft and hard 
polydomains has not previously been appreciated, but 
very recent experiments confirm that it is correct [Io| . 
The recognition of softness in some polydomain systems 
makes the fabrication of useful LCE soft actuators more 
likely since polydomains are much easier to synthesize 
and arc not limited to thin film geometries. Our results 
suggest that a SmC* polydomain cross-linked in a SmA 
monodomain state would be a good choice for low field 
electrical actuation. 

LCE's are very analogous to martensitic metals since 
both systems exhibit symmetry-breaking transitions cou- 
pled to deformations. In the martensitic case the sym- 
metries that are broken are discrete whereas in LCE's 
they are continuous. However, drawing analogies be- 
tween LCE polydomains and martensite polycrystals is 
quite subtle. The soft polydomain LCE's with homo- 
geneous high temperature cross-linking states are analo- 
gous to single crystal martensite systems. However, sin- 
gle martensite crystals arc difficult to prepare because, 
even if they are prepared in the high symmetry state, 
they are not isotopic so the crystal can form with dif- 
ferent orientations at different points in space. This 
is in marked contrast to isotropic genesis LCE poly- 
domains where, because the cross-linking state is com- 
pletely isotropic, cross-linking in a spatially homogeneous 
state is trivial. There is no satisfactory martensite ana- 
log of cross- linking in the low symmetry polydomain state 
since martensite poly crystals are formed in a poly crys- 
talline high symmetry state, so they access a stress-free 
high symmetry state, whereas the low symmetry cross- 
linked elastomers cannot. There is one, albeit rather con- 
trived, LCE system, not analyzed here, that is directly 
analogous to a martensite polycrystal — an elastomer 
cross-linked in the SmA polydomain state then cooled to 
a SmC polydomain state. Such an elastomer would, on 
heating, return to a stress-free SmA polydomain state, 
and has local soft modes in the SmC state generated by 
the symmetry breaking SmA-SmC transition. 
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